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Abstract. A Quantum Computer is a new type of computer which can 
solve problems such as factoring and database search very efficiently. 
The usefulness of a quantum computer is limited by the effect of two 
different types of errors, decoherence and inaccuracies. In this paper we 
show the results of simulations of a quantum computer which consider 
both decoherence and inaccuracies. We simulate circuits which factor 
the numbers 15, 21, 35, and 57 as well as circuits which use database 
search to solve the circuit satisfaction problem. Our simulations show 
that the error rate per gate is on the order of 10~ 6 for a trapped ion 
quantum computer whose noise is kept below 7r/4096 per gate and with 
a decoherence rate of 10 -6 . This is an important bound because previous 
studies have shown that a quantum computer can factor more efficiently 
than a classical computer if the error rate is of order 1CP 6 . 
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1 Introduction 

A quantum co mputer consists of atomic particles which obey the laws of quan- 
tum mechanics TuHo95 1 1 Lloy9Ej| . The complexity of a quantum system is expo- 



nential with respect to the number of particles. Performing computation using 
these quantum particles results in an exponential amount of calculation in a 



polynomial amount of space and time |Feyn85||Deut85|. This quantum paral- 



lelism is only applicable in a limited domain. Prime factorization is one such 



problem which can make effective use of quantum parallelism! Shor94|. This is 



an important problem because the security of the RSA public-key cryptosystem 



relies on the fact that prime factorization is computationally difficult [ PJSA78 1 . 

Errors limit the effectiveness of any physical realization of a quantum com- 
puter. A quantum computer is subject to two different types of errors, decoher- 
ence and inaccuracies. Decoherence occurs when a quantum computer interacts 
with the environment. This interaction destroys the quantum parallelism by 
turning a quantum calculation into a classical one. The other type of error, in- 
accuracies in the implementation of gate operations, accumulates over time and 
destroys the results of the calculation. 

In this paper we show results of simulations of a quantum computer which 
is subject to both decoherence and inaccuracies. These simulations assume the 



trapped ion model of a quantum computer proposed by Cirac and Zoller[CiZo9£ ] . 



We study Shor's factorizati on algor ithm by simulating circuits which factor the 



numbers 15, 21, 35, and 57| Shor94 |. We also simulate Grover's databa se searc h 



algorithm with a circuit which solves the circuit satisfaction problem |Groy_9£]. 
The rest of this section gives a brief overview of quantum computers. 

1.1 Qubits and Quantum Superposition 

The basic unit of storage in a Quantum Computer is the qubit. A qubit is like 
a classical bit in that it can be in two states, zero or one. The qubit differs 
from the classical bit in that, because of the properties of quantum mechanics, 
it can be in both these states simultaneously! FeLS65( | . A qubit which contains 



both the zero and one values is said to be in the superposition of the zero 
and one states. The superposition state persists until we perform an external 
measurement. This measurement operation forces the state to one of the two 
values. Because the measurement determines without doubt the value of the 
qubit, we must describe the possible states which exist before the measurement 
in terms of their probability of occurrence. These qubit probabilities must always 
sum to one because they represent all possible values for the qubit. 

The quantum simulator represents the qubits of the computer using a com- 
plex vector space. Each state in the vector represents one of the possible values 
for the qubits. The bit values of a state are encoded as the index of that state 
in the vector. The simulator represents each encoded bit string with a non zero 
amplitude in the state vector. The probability of each state is defined as the 



square of this complex amplitude [ FeL S 6 5 1 . For a register with M qubits, the 
simulator uses a vector space of dimension 2 M . 

1.2 Quantum Transformations and Logic Gates 

A quantum computation is a sequence of transformations performed on the 



qubits contained in quantum registers! Fcyn85 [BaBe95|. A transformation takes 



an input quantum state and produces a modified output quantum state. Typ- 
ically we define transformations at the gate level, i.e. transformations which 
perform logic functions. The simulator performs each transformation by multi- 
plying the 2 dimensional vector by a 2 x 2 M transformation matrix. 

The basic gate used in quantum computation is the controlled-not, i.e. exclu- 
sive or gate. The controlled-not gate is a two bit operation between a control bit 
and a resultant bit. The operation of the gate leaves the control bit unchanged, 
but conditionally flips the resultant bit based on the value of the control bit. 

1.3 The Ion Trap Quantum Computer 

The ion trap quantum computer as proposed by Cirac and Zoller is one of the 
most promising schemes for the experimental realization of a quantum com- 
puter JCiZo95|]. Sever al experiments have demonstrated simple quantum gates 



MoMe95] [ WiMM96 ] . Laser pulses directed at the ions in the trap cause trans- 



formations to their internal state. A common phonon vibration mode is used to 



communicate between the ions in the trap. A controlled-not gate is a sequence 
of laser pulses. We use the ion trap quantum computer as the model for our 
quantum simulator. 



Qubits in the Ion Trap Quantum Computer. Qubits are represented using 
the internal energy states of the ions in the trap. The ion trap represents a logic 
zero with the ground state of an ion, and a logic one with a higher energy 
state. The ion trap quantum computer also requires a third state which it uses 
to implement the controlled-not gate. In this paper we use a simplified model 
which, instead of using a third state for each qubit, uses a single third state which 
is shared amongst all the qubits. This simplified model reduces t he simula tion 
complexity exponentially without an appreciable loss of accuracy [ ObDe97a | . 



Transformations in the Ion Trap. An operation in the ion trap quantum 
computer is a sequence of laser pulses. Each laser pulse is defined by one of the 
transformation matrices shown in equation |l|. 9 corresponds to the duration of 
the laser pulse and (j> corresponds to the phase. A two bit controlled-not gate is 
a sequence of five laser pulses, two V and three U transformations. A single bit 
not gate can be implemented with three laser pulses and the three bit controllcd- 
controlled-not gate requires seven laser pulses. 
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2 Simulating a Quantum Computer 

Our quantum computer simulator simulates circuits at the gate level. The simu- 
lator implements one, two and three bit controlled-not gates as well as rotation 
gates. The simulator implements each gate as a sequence of laser pulses, and 
represents the entire vector space throughout the simulation. 



2.1 Operational Errors and Decoherence 

The simulator models inaccuracies by adding a small deviation to the two angles 
of rotation 9 and <j). Each operational error angle is drawn from a gaussian 
distribution with a parametrized mean (/i) and standard deviation (a). Errors 
with non zero ^i, called mean errors, correspond to systematic calibration errors, 
and errors with non zero <r, referred to as standard deviation errors, correspond 
to noise in the laser apparatus. 

Because the phonon mode is coupled to all the qubits in the computer, it is 



the largest source of decoherence [MoMe95 . For this reason we only model the 



phonon decoherence and not the decoherence of the individual qubits. 



We model the decoherence of the phonon mode by performing an additional 
operation after each laser pulse. Equation || shows this transformation which 
has the effect of decaying the amplitude of the states in the phonon state. This 



decay transformation is based on the quantum jump approach | Carm93 . The 

(2) 



decay parameter (dec) remains constant throughout the entire simulation. 

|V > |o > p \%i> > |o > P 



\iP > |1 > p e~ dec / 2 \4> > |1 > p 

This method of modeling decoherence implicitly models spontaneous emission. 
Because the state is never renormalized, the total norm at each step represents 
the probability that the calculation survives up to that point without a sponta- 
neous emission occurring. An alternative method for modeling decoherence is to 
renormalize the state at each step and then, based on a probability of emission, 
cause emissions at different points in the calculation. This method has the dis- 
advantage that, because we cause emissions at random points in the calculation, 
we must run multiple simulations each with different initial random seeds to 
average out any bias caused by the random number generator. We have shown 
however that both methods for modeling decoherence give essentially the same 



results [ObDe97a] 



Because the simulator applies the decoherence transformation once per laser 
pulse, the parameter dec has units of (decoherence/laser pulse). To convert these 
units to decoherence per unit time we must consider the switching time of the 
laser. The dec parameter is simply the switching speed divided by the decoher- 
ence time. Recent experiments show switching speeds of 20kHz for a controlled- 
not gate, i.e. four ir pulses, and a decoherence rate of a few kHz[ MoMe95|] . This 
corresponds to a decoherence parameter value between 10~ 2 and 10 



-3 



2.2 Quantum Circuits 

Much of the current interest in quantum computation is due to the discovery 
of an efficient algorithm by Peter Shor to factor numbers [ [5hor94 . By putting 



the qubit register A in the superposition of all values and calculating the func- 
tion f(A) = X A modN, a quantum computer calculates all the values of f(A) 
simultaneously. Where N is the number to be factored, and X is a randomly 
selected number which is relatively prime to N. The quantum factoring circuit 
also contains operations to create the superposition state at the beginning of the 
circuit, and extract the period at the end of the circuit. The circuit to calculate 



f(A) can be performed in 0(L 3 ) time using repeated squaring} Desp96 |, i.e. a 
sequence of multiplications performed modulo N. 

Grover's database search algorithm searches for a key, from a set of match- 



ing keys, in an unsorted database! Grov96 . The keys are defined by a function 



which can be evaluated in unit time. After evaluating this function a diffusion 
transformation is performed which amplifies the probability in the states with 
matching keys. Grover shows that after performing 0(y/(N)) evaluation steps 
and diffusion transformations the probability of measuring a matching key is 
greater than 1/2. 



3 Simulation Results 



In this section we study how decoherence and operational errors degrade the 
fidelity of the factoring and database search algorithms. 

The fidelity as defined by fidelity = \\ < <yj|V> > || 2 , measures how close a 
state with error in it is to the correct result. The fidelity is defined as the inner 
product between the simulation with errors (jp) and the correct result (</?). 

Table 1 shows all the benchmarks used in our simulation studies. To show 
the complexity of simulating these benchmarks we show the simulation time for 
each. This time assumes a single 300 MHz processor and the simulations include 
only operational error. All the factoring benchmarks were run using a parallel 
version of the simulator! ObDc9S|. We have run simulations on as many as 256 
processors, and the simulator achieves near linear speedup. 



Table 1. Benchmarks used in simulation studies 





Number of 


Number of 




Simulation 


Benchmark 


qubits 


laser pulses 


Description 


time (sees) 








Circuit SAT using the Grover 




grover 


13 


1,838 


database search algorithm 


10 








One modulo multiply step 




mult 


16 


8,854 


from the factor- 15 problem 


282 








Factor- 15 problem using 




factorl5 


18 


70,793 


3 qubits for A 


10,465 








Factor-21 problem using 




factor21 


24 


69,884 


6 qubits for A 


272,276 








Factor-35 problem using 




factor35 


27 


99,387 


6 qubits for A 


3,083,520 








Factor-57 problem using 




factor57 


27 


97,939 


6 qubits for A 


3,067,853 



Because of round off error, in the factoring algorithm, the choice of the num- 
ber of bits to use in the A register affects the probability seen after performing 
the FFT[Shor94|. Shor suggest using 2L + 1 bits for an L bit factorization, but 
for the numbers used in our studies we can use less. For the factorl5 problem the 
period is a power of two, i.e. four, and therefore there is no round off error. For 
the numbers 21, 35 and 57, the probability given by the FFT does not increase 
for more than six bits. Also we are mainly concerned with observing the fidelity 
for these circuits, and the fidelity is always calculated before the FFT. 



3.1 Operational Errors 

In this section we consider operational errors without any decoherence. We first 
investigate the significance of errors in the angle <j), by varying the amount of 



error in the angles 9 and <fi separately. In all other simulations we vary the angles 
9 and <f> together. 

To average the random bias out of simulations with non zero a we run mul- 
tiple simulations each with different initial random seeds. To get a single sim- 
ulation point we run at least four simulations, and in many cases we run more 
to establ ish uppe r and lower confidence intervals for the average of the ending 
fidelities | H1M08O | . For simulations which include errors in the angle 9 we run 
enough simulations so that the average of the fidelities is within 0.02 of the ac- 
tual mean of the distribution with an upper and lower confidence of 95%. There 
is more variability in the fidelity for simulations which consider only <f> errors, so 
for these we obtain the 95% confidence intervals to within 0.03 of the mean. 



Significance of Errors in the Angle <fi. Fig. 1 shows how operational errors 
degrade the fidelity for the factor 15 benchmark. In Fig. 1(a) we vary the mean 
and standard deviation and introduce both 9 and (f> errors. In Fig. 1(b) we only 
introduce 9 errors. Comparing the two graphs shows that the combination of 9 
and <f> errors produces a lower fidelity than 9 errors alone. 




Fig. 1. Inaccuracies for factorl5 (a) 9 and <f> errors (b) 9 errors (c) (j> errors (d) 9 and 
<j> errors, 9 errors, 4> errors 



Fig. 1(c) shows the effect of only errors. As the graph shows noise degrades 
the fidelity, but adding a constant amount of error has no effect. Fixed magnitude 



<j> errors have no effect because the laser transformations are always performed 
in pairs, and an error in the second transformation cancels an error in the first 
transformation. In Fig. 1(d) we compare the effect of 9 and <f> errors alone and 
their combined effect. The highest degradation occurs when considering both 9 
and 4> errors, and 9 errors produce a more significant effect than <f> errors. 



Operational Errors for the Grover Benchmark. Fig. 2 shows mean and 
standard deviation operational errors for the grover benchmark. The figures show 
the probability of finding a correct key for twelve iterations of the algorithm. 
When running the database search it is important to stop after the correct iter- 
ation, because the probability decreases if we run too many iterations BoBr96|. 
For our case the highest probability occurs after the eighth iteration. 



(a) Grover, (I = 0, No Decoherence 
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Fig. 2. Inaccuracies for the grover benchmark, (a) a Errors, (b) fi errors 



Fig. 2(a) shows that the peak probability is above 0.5 for a errors as great 
as 7r/128, and for an error rate of ir/64 the peak probability is 0.2. Also for 
it errors the peak probability occurs after the same iteration for all errors less 
than 7T/T28. However // errors, as shown in Fig. 2(b), shift the peak so that it 
occurs at an earlier iteration. The peaks are higher for \x errors than they are 
for a errors with the same level of error. However, since we can only perform a 
single measurement, the shift in the peak values causes a further reduction in 
the probability if we measure after the eighth iteration. 



Fidelity at Intermediate Points for the Factoring Benchmarks. Fig. 3 
shows the fidelity at intermediate points in the calculation for the factor21, and 
factor57 benchmarks. 

Standard deviation errors produce a more significant effect than /i errors of 
the same magnitude. This is due to a cancellation effect for mean errors which 
is very similar to the cancellation effect exhibited by 4> errors pbDc97a| . This 
cancellation effect occurs because of the nature of reversible computation and 
causes the fidelity to go through periods where it increases. Because we perform 
operations in pairs, where the two operations are inverses of each other, an 
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(b) Factor2 1 , (l = 0, No Decoherence 
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(d) Factor57, u. = 0, No Decoherence 
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Fig. 3. Fidelity as a function of the number of gates for factor21 and factor57 



error in one operation may reverse the error from an earlier operation. Standard 
deviation errors also exhibit this effect but it is reduced because we average the 
results from multiple simulations. 



3.2 Decoherence Errors 

Fig. 4 shows how the fidelity of a computation decreases over time in the presence 
of decoherence. Fig. 4 shows the fidelity at intermediate points in the calculation 
for the four factoring benchmarks. For small amounts of decoherence, i.e. 10~ 6 or 
less, the fidelity is not adversely affected. A decoherence rate of 10~ 4 results in a 
steady decrease in the fidelity over the course of the computation, and for rates 
even higher the fidelity drops off very quickly. As the figure shows, decoherence 
has a similar effect on all of the benchmarks. 



3.3 The Correlation Between Decoherence and Operational Errors 

Both decoherence and operational error cause a degradation of the fidelity in a 
quantum computation. Decoherence degrades the fidelity through the decay of 
the phonon state, and operational errors result in the accumulation of amplitude 
in unwanted states. The combined effect of these two factors is a degradation 
which is worse than either factor considered alone. We can represent the com- 
bined effect as: F dec . op = F dec • F op + S7(F dec , F op ) Where F dec and F op arc the 



(a) Factorl5, Decoherence, No Operational Errors 




(b) Factor21, Decoherence, No Operational Errors 
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Fig. 4. Fidelity as a function of the number of gates for decoherence in the factorl5, 
factor21, factor35 and factor57 benchmarks 



fidelities of simulations for decoherence and operational error considered sepa- 
rately, and Sl(Fdec, F op ) is the correlation between the two types of error. As 
Table 2 shows the correlation is very low. We calculated the correlation by run- 
ning simulations which considered decoherence and operational error together. 
For all the benchmarks the maximum correlation is at most 1.14 x 10~ 2 . This 
result means that we can simulate decoherence and operational errors separately, 
and combine the results to obtain their collective effect on a calculation. 



Table 2. Correlation (H) between decoherence and operational errors 



Benchmark and Simulation Model 


Maximum Q 


Average Q 


mult, ^=0, a = tt/1024 - tt/64 


5.76 x 10" 5 


3.63 x 10~ 6 


mult, a = 0, fi = tt/1024 - tt/64 


9.26 x 10" 3 


5.51 x 10~ 4 


factorl5, a = tt/1024, /x = 


4.15 x 10" 3 


3.96 x 10~ 4 


factorl5, a = 0, fi = tt/1024: 


1.14 x 10~ 2 


8.76 x 10~ 4 


grover, n=0, a = tt/1024 - tt/128 


1.78 x 10" 3 


1.02 x 10 -4 


grover, cr=0, fi = tt/1024 - tt/128 


2.67 x 10" 3 


2.36 x 10~ 4 



The complexity of simulating decoherence alone is much lower than the com- 
plexity of simulating operational errors. This is because the decay transformation 
does not have any off diagonal terms, and therefore it does not introduce any 
new error states. The simulator only needs to represent enough states to repre- 
sent the superposition state. Instead of using the index of a state to represent 
its bit value, the simulator now keeps an extra field for each state which holds 
the current value of the bit string for that state. 

Using this new method for modeling decoherence, the simulator only needs 
to allocate 0(2 L ) states to simulate the factorization of an L bit number. To 
represent all the qubits for this problem, the simulator would need to allocate 
0(2 4L ) states. This reduces the memory requirements and simulation complexity 
by a factor of 0(2 3L ). For example the simulation of the factorl5 circuit requires 
only 1/4096 the amount of time as before. 



3.4 The Error Rate per Gate 

As our simulation results show, a modest amount of error destroys even a rel- 
atively small calculation. If quantum computers are to be useful, we must be 
able to perform calculations which are even larger than the ones considered 
here. These larger calculations will therefore require the use of quantum error 



correcting codes [3tea96 . Several recent studies have shown that, by using fault 
tolerant techniques, if the error of an individual gate is low enough we can per- 
form a useful quantum calculation of indefinite length [ Pres96| | KnLZ96 ] . This 



accuracy threshold is expressed in terms of the probability of error per gate. We 
can use the results of our simulation studies to show how this error probability 
relates to decoherence and inaccuracies. 

Table 3 shows the error rate per gate considering decoherence and operational 
errors for the factor57 benchmark. Error rates for the other factoring benchmarks 
as well as error rates for the combination of decoherence and operational errors 



are given in | ObDc97b |. The error rate is calculated as (1 - Fidelity) /Number of 
Gates. To get the error rate for a particular amount of error, we calculate the 
error rate after every tenth gate in a computation and take the average. A gate 
is either a one, two or three bit controlled-not gate or a single bit rotation. It 
takes five laser pulses on average to implement each gate. 

To perform a computation of arbitrary length the error rate must be about 
10~ 5 for one and two bit gates and 10 -3 for three bit gates [Pres96]. An error 
rate of 10~ 5 corresponds to operational errors with a = 7r/2048 and decoherence 
of 10~ 5 . Table 3 shows that we can tolerate an even higher level of constant 
magnitude errors. The error rate is 10 -5 for operational errors with /i = 7r/1024. 

To perform a quantum factorization which is more efficient than a classi- 
cal one the error threshold is even tighter, roughly 10~ 6 . This corresponds to 
operational errors with a < 7r/4096 and a decoherence rate of 10~ 6 . 

Using the error rate of our factoring circuits to predict the error rate for 
error correction circuits assumes that these two types of circuits behave in a 
similar manner. These circuits are similar because they are both built from the 
same types of elementary gates, i.e. controlled-not gates. Also the error rate, 



Table 3. Average error rate per gate for the factor57 benchmark 



Decoherence 


Operational Errors, fi = 


Operational Errors, a — 


dec 


Error Rate 


a 


Error Rate 


A* 


Error Rate 


ltr 7 


9.1 x 1(T 8 


tt/8192 


6.6 x 1(T 7 


tt/4096 


7.9 x 10~ 7 


ltr 6 


9.1 x 10~ 7 


71-/4096 


2.6 x 10~ 6 


tt/2048 


3.1 x 10~ 6 


10- 5 


8.8 x 1(T 6 


tt/2048 


1.0 x 10~ 5 


tt/1024 


1.2 x 10 -6 


ltr 4 


6.5 x 1(T 5 


tt/1024 


3.6 x 10~ 5 


tt/512 


4.4 x 10~ 5 



for a given amount of error, is very similar for all the factoring circuits that 
we considered. Lastly effects that we have seen such as error cancellation are 
a by-product of the fact that quantum circuits and gates are implemented in a 
reversible fashion. Because of the nature quantum mechanics quantum circuits 
will always be implemented in this way. 

4 Conclusion 

Quantum computation is a new type of computation which can achieve expo- 
nential parallelism. The feasibility of a quantum computer is threatened by two 
types of errors, decoherence and inaccuracies. In this paper we performed simu- 
lations of a quantum computer running Shor's factoring algorithm and Grover's 
database search algorithm to access the feasibility of implementing a quantum 
computer. 

Our simulations show that random inaccuracies (noise) are more significant 
than fixed magnitude inaccuracies for the ion trap quantum computer. Also 
errors in the duration of the laser pulse are more significant than errors in the 
phase of the laser. For the problems considered in this paper we show that 
noise or constant magnitude inaccuracies of magnitude 7r/512 or greater cause 
a significant impact on the fidelity of the calculation. 

Our simulations also show that a quantum computation can tolerate a deco- 
herence rate as high as 10~ 5 . For the ion trap computer this corresponds roughly 
to a decoherence lifetime of 50 milliseconds and a switching speed of 500 nanosec- 
onds. We also show that inaccuracies and decoherence are uncorrelated and can 
be simulated separately. 

Our simulations relate the physical quantities of inaccuracies and decoherence 
to the probability of error per gate. An error rate per gate on the order of 
10~ 6 corresponds to inaccuracies of less than 7r/4096 per laser operation and a 
decoherence rate of 10 -6 . A quantum computer with this error rate, if it uses 
quantum error correcting codes, could factor a number more efficiently than a 
classical computer. This assumes that the quantum circuits used to implement 
factoring with error correction codes behave in the same manner as the factoring 
circuits used in this paper. 
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